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Abstract 

Inclusive pion production at high transverse momenta (px ^ 2 GeV/c) in nucleus- 
nucleus (A+A) collisions at CERN SPS {^/s NN « 20 GeV) is revisited and system- 
atically compared to all existing proton-proton data in the same range of center-of- 
mass energies. The ratio of A+A to p+p pion cross-sections (nuclear modification 
factor) for central Pb+Pb, Pb+Au and S+Au reactions does not show a strong 
enhancement as a function of pr as previously found, but is consistent with scaling 
with the number of nucleon-nucleon (NN) collisions. Neutral pion yields per NN 
collision in head-on Pb+Pb reactions are suppressed, whereas peripheral yields are 
enhanced. These results together indicate that some amount of "jet quenching" may 
already be present in central heavy-ion reactions at Js w 20 GeV. 



PACS: 25.75.-q, 12.38.Mh, 13.85.-t, 13.87.Ce 



1 Introduction 



Lattice calculations of bulk Quantum Chromodynamics (QCD) in thermal 
equilibrium [1] predict the transition of hadronic matter to a deconfmed and 
chirally symmetric system of quarks and gluons above energy densities of the 
order e cr j t ~ 0.7 ± 0.3 GeV/fm 3 . The formation and study under laboratory 
conditions of this "Quark-Gluon Plasma" (QGP) phase is one of the highest 
priorities in high-energy nuclear physics in the present day. Several experimen- 
tal results from the CERN Super Proton Synchrotron (SPS) relativistic heavy- 
ion programme collected during the 1990's in fixed-target experiments with 
center-of-mass energies y/s NN ~ 20 GeV have been interpreted, not without 
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controversy, in terms of QGP formation [2] . Indeed, although several observa- 
tions in central Pb+Pb collisions [3] are consistent with expected QGP signals, 
e.g. the "anomalous" suppression of charmonium states due to Debye screening 
of the color potential in the plasma [4] (see, however, also [5]), other expected 
signatures such as "jet quenching" due to parton energy loss in the dense de- 
confined medium [6] seem to be significatively absent from the data. In fact, 
high pt pion production in central A+A at CERN-SPS was found not to be 
suppressed but enhanced compared to production in free space [7,8,9]. Such a 
"Cronin enhancement", observed earlier in p+A [10,11,12] and a + a [13] colli- 
sions, is usually interpreted in terms of multiple initial-state parton scatterings 
which result in a broadening of the final Pt spectra [14]. In contrast, high-p^ 
hadro-production in central Au+Au collisions at y/s NN = 130 GeV [15,16] and 
200 GeV [17,18,19,20] at the BNL Relativistic Heavy Ion Collider (RHIC) has 
been found to be strongly suppressed (by up to a factor of 4-5) compared to 
p+p collisions measured at the same yj s NN [21,18]. The observed suppression 
at RHIC is even more significant considering the fact that the "default" pro- 
duction at high pt in the "cold" nuclear environment of deuteron-nucleus re- 
actions at collider energies is also Cronin enhanced [22,23,24,20]. These results 
clearly indicate that final-state effects are responsible for the high p T deficit 
observed in central Au+Au. The difference between the suppressed RHIC and 
enhanced SPS hadro-production at large pt implies that there must exist an 
intermediate value of collision energies in nucleus-nucleus reactions at which 
final-state quenching starts to dominate over initial-state Pt broadening. The 
search for the onset of high p? suppression is the main motivation behind the 
dedicated v /s^ = 62.4 GeV Au+Au run carried out at RHIC in April 2004. 

Theoretical studies [25] of parton propagation in a dense medium show that 
the induced parton energy loss is proportional to the initial gluon density 
(dN 9 /dy) in the system. In this context, the absence of suppression at CERN- 
SPS is surprising considering that experimental estimates of the initial energy 
density, based on the Bjorken prescription [26] for a boost-invariant longitu- 
dinally expanding plasma, at a canonical thermalisation time of tq = 1 fm/c, 
are on the order of e|f (r ) « 3 GeV/fm 3 [27,28] and e«f IC (r ) « 5 GeV/fm 3 
[29]. Such energy densities, both well above the critical value of ~0.7 GeV/fm 3 , 
would correspond to a ratio of parton rapidity densities (dN/ dy oc p oc e 3 / 4 for 
an ideal gas of quarks and gluons) of ~0.68 between RHIC and SPS at r = 1 
fm/c. Since -R™ c ~ 0.2, one would accordingly expect suppression factors of 
the order R s /X ~ 0-3 assuming the same Cronin pr broadening at RHIC and 
SPS (see discussion later). This is in violent contradiction with the R s /X >>1 
usually quoted in the literature. The usual explanations for the reported ab- 
sence of suppression at SPS involve short QGP lifetime [8,9], domination of 
multiple soft collisions (Cronin effect) over hard scatterings [30], and modest 
amount of parton rescatterings [31]. 

In this letter we explore an alternative interpretation based on a thorough re- 
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analysis of the p+p — > 7r°+X baseline spectrum used to determine the A+A 
nuclear modification factor at y/s = 17.3 GeV in [7,9]. It turns out that the 
absence of a concurrent proton-proton measurement at the same y/s, and the 
use of inexact baseline references extrapolated from higher collision energies, 
result in apparent strong Cronin enhancements which are not actually sup- 
ported by the data. The fact that the p+p reference spectra for SPS used 
so far are not well under control can be already realized by inspecting the 
original work [7] which shows two results difficult to reconcile at first sight: a 
suppressed n° production as given by the ratio of central to peripheral Pb+Pb 
collisions (R cp ~ 0.6), and a strongly enhanced central Pb+Pb over p+p ratio 
(albeit with large systematic uncertainties), Raa ~ 8 at the largest p?- In 
this paper, the whole set of available high p? data from the SPS heavy-ion 
experiments: ir° and 7r =l= at y/s NN = 17.3 GeV from Pb+Pb (WA98) [7] and 
Pb+Au (CERES/NA45) [32] respectively, and tt° from S+Au at y/I^ = 19.4 
GeV (WA80) [33], will be reexamined and carefully compared to p+p spec- 
tra constructed from a ^/s-dependent global fit of most of the available pion 
differential cross-sections in the range y/s ~ 7 - 63 GeV [34]. Using this new 
reference, it will be shown that high-py hadropro duct ion at y/s NN ~ 20 GeV is 
not enhanced in central nucleus- nucleus reactions but, within errors, is consis- 
tent instead with scaling with the number of NN collisions. Furthermore, the 
fact that high-pT vr° production in the the top 1% central lead-lead collisions 
appears to be suppressed by a factor of ~1.6 ± 0.6 compared to the p+p 
reference and that the peripheral yields are Cronin enhanced, points to some 
mechanism of final-state suppression at work in central Pb+Pb at CERN-SPS 
energies. The revised nuclear modification factors will be compared to the pre- 
dictions of a pQCD-based model of parton energy loss, and the (partonic or 
hadronic) nature of the dissipative medium will be discussed. 



2 High p T pion production in p+p collisions at y/s 20 GeV 

Particle production at high px {pr ^ 2 GeV/c) in hadronic collisions re- 
sults from incoherent parton-parton scatterings with large Q 2 . In the ab- 
sence of initial and final state interactions, independent scattering and pQCD 
factorization 1 imply that inclusive A+B cross-sections for hard processes 
should scale simply as A ■ B times the corresponding p+p cross-sections: 
E da^ard/ d 3 P = A • B • E da^ard/d^P- Usually heavy-ion experiments mea- 
sure invariant yields for a given centrality bin and, thus, the corresponding 
"scaling law" reads EdN^ d /d 3 p = (T AB (b)} ■ EdaZ r d/d 3 P, where T AB (b) 

1 Incoming quarks and gluons undergoing hard scattering are considered "free" in 
a collinear factorized approach, i.e. the density of partons in a nucleus with atomic 
number A is considered to be equivalent to the superposition of A independent 
nucleons, or f a /A( x >Q 2 ) = A f a / p (x,Q 2 ) in terms of parton distribution functions. 
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is the Glauber nuclear overlap function at impact parameter b [35]. Since 
the number of inelastic nucleon-nucleon collisions at b is proportional to Tab 
(N coll (b) = T AB {b) ■ a™ e \ with a l p n p el = 32 mb at y/s « 20 GeV), one usually 
quantifies medium effects at high pt via the nuclear modification factor 

t? f \ —— d 2 N AB /dydp T 
ab (Pt) - {Ncou{b)) x d 2 Npp/dydpT ' I ) 



which measures the deviation of A+B at impact parameter b from an in- 
coherent superposition of nucleon-nucleon collisions. Unfortunately, on the 
experimental side no p+p — > n°+X spectra measured at high p? close to mid- 
rapidity exists at the same collision energy 2 of the SPS Pb-induced nuclear 
collisions (lfj & = 158 AGeV corresponding to y/s = 17.3 GeV). On the the- 
oretical side, the single inclusive particle spectrum at large pt in high-energy 
p+p collisions can be in principle calculated within the framework of collinear 
factorization. However, below y/s ps 60 GeV as noted already in [39], the 
cross-sections at pr < 5 GeV/c in p+p collisions are underpredicted by stan- 
dard pQCD calculations, and additional non-perturbative effects (e.g. intrinsic 
kr) must be introduced to bring parton model analysis into agreement with 
data. Unfortunately, those effects cannot so far be introduced in a model- 
independent way and different pQCD calculations [8,30,38] effectively yield 
different final pion spectra for p+p collisions around y/s = 20 GeV. 

In the absence of a concurrent experimental measurement and lacking a fully 
reliable theoretical calculation, two approaches have been followed to construct 
a p+p reference for n° production at SPS energies. First, the WA98 collabo- 
ration [7] has employed a semi-empirical modified power-law form A [p / (pr + 
Po)] n (originally proposed by Hagedorn [40]) tuned to reproduce the pt spec- 
tra measured at higher y/s, plus an xt scaling prescription [41] to account for 
the collision energy dependence of the cross section. Second, Wang&Wang [9] 
have adopted a more complex power-law ansatz for the pt spectrum which de- 
scribes the charged pion data at y/s = 19.4 GeV [11], combined with a pQCD 
parton model calculation to scale the cross-section down to y/s = 17.3 GeV. 
Both parametrizations have been tuned to reproduce a subset of the existing 
p+p — > 7r+X data at y/s ~ 20 GeV, but no true global analysis has been 
carried out to fully compare the parametrizations to all the existing results 
in this energy regime. Figure 1 shows a comparison of parametrizations [7] 
and [9] to the whole set of experimental results on high p? pion production 
at 9 cm ~ 90° in the y/s = 16.9 - 19.4 GeV range. The first thing worth to 
notice is the relatively large disparity among the experimental data obtained 

2 Ref. [36] provides a p+p — > tt°+X measurement at pi a t = 158.9 GeV/c very close 
to WA98 Pb+Pb beam energy. However, these data are inconsistent with the rest 
of the tabulated results (the authors themselves have discarded this sample from 
subsequent analysis [37]) and so have not been considered in this work. 
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at the same energy (y/s = 19.4 GeV), especially at the largest values. This 
fact highlights the importance of concurrently measuring A+A and baseline 
p+p differential cross-sections at the same center-of-mass energy and with the 
same setup in present and future heavy-ion experiments. Having said that, 
the WA98 reference fit [7] reproduces the available data only at p T « 2 GeV/c 
whereas it systematically overpredicts the yields below this px value and signif- 
icantly underpredicts them above it. Wang&Wang parametrization [9], repro- 
duces better the shape of the pt spectra but it systematically underpredicts 
most of the ir° yields by ~50% in the p T xs 1.5 - 4.0 GeV/c range. Both 
data fits, however, seem to follow rather closely the (n + + n~)/2 data 3 from 
Antreasyan et al. [11] which is actually used as the basic set for constraining 
the fit parameters at y/s = 19.4 GeV in [9]. Forcing the parametrization [9] to 
fit the 7r =l= data of [11] above 4 GeV/c (where the scarce A+A data have large 
statistical errors anyway) without any other constraint from existing n° results 
in the same range of collision energies seems to be the cause of the limited 
agreement of this parametrization with the proton-proton data in the range 
Pt ~ 1.5 - 4.0 GeV/c where the heavy-ion measurements are available. In 
the case of the parametrized 7r° reference of WA98, the disagreement data-fit 
seems to come from the ansatz used to take into account the y/s dependence 
of the cross section since the original parametrization reproduces well the ISR 
7t° spectra at higher energies [7]. 

At variance with these two works, Blattnig et. al [34] derive a parametriza- 
tion of the invariant differential cross-section for inclusive ir° production in 
proton-proton collisions based on a global analysis of most of the available 
data within y/s ^7-63 GeV. This is a purely empirical 11-parameter func- 
tional form tuned to provide a reasonable description of the full spectral 
shape, angular distribution and total cross-section of the measured pions in 
this range of center-of-mass energies. Fig. 2 shows the level of agreement of 
Blattnig fit to the available pion data in p+p collisions for beam energies in 
the range of CERN-SPS heavy-ion experiments. The agreement between data 
and fit is more satisfactory, both in shape and magnitude, than the two pre- 
vious parametrizations especially within px ~ 1.5 - 3.5 GeV/c. It describes 
rather well, in particular, the most recent (and precise) data sets from Fer- 
milab E704 experiment [43] which were not actually considered in the fitting 
analysis of [34]. It is worth noticing, however, that at pt > 3.5 GeV/c al- 
though the parametrization reproduces the 7r° data sets of ref. [36], it seems 
to be ~50% above the n° results of refs. [42,43] and the averaged 7r ± of [11]. 
Nonetheless, given the relatively poor agreement between p+p experimental 



3 Pion production in this kinematical range of parton fractional momenta in the 
proton, x > 0.2, is dominated by valence u (2 quarks) and d (1 quark) fragmentation 
which, due to their respective quark content, results in the o(t[ + ) < a(ir°) < cf(it~) 
ordering of the pion cross-sections. One would, thus, indeed expect (ir + + vr~)/2 to 
provide a good approximation of the 7r° yields in p+p collisions. 
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Fig. 1. Relative differences between the single inclusive pion spectra measured in 
p+p collisions at y/s = 16.9 - 19.4 GeV [11,36,42,43] and the p+p -► ir°+X 
parametrizations proposed by the WA98 collaboration [7] (upper figure) and 
Wang&Wang [9] (lower figure) at the corresponding ^/s. The shaded band rep- 
resents the 20% overall uncertainty originally assigned to the WA98 parametriza- 
tion. Wang&Wang only provides the fit parameters for 2 fixed values of = 17.3, 
19.4 GeV. 

measurements themselves above pr ~ 3.5 GeV/c, and especially given that the 
highest pr vt° measured in A+A reactions are at ~4 GeV/c (with statistical 
errors which are larger than the p+p reference uncertainty), we consider this 
fit to provide a much more accurate representation of the proton-proton 7r° 
baseline production than the two previously used parametrizations. Hereafter, 
we will thus use the reference of Blattnig et. al [34] as a reasonable estimate 
of the p+p^ 7T°+X spectra for ^Js = 17-20 GeV in the range p T — 1 - 
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4 GeV/c, - with an assigned overall uncertainty of ±25% (shaded band in 
Fig. 2) to cover most of the existing measurements as a benchmark for the 
study of high p T pion production in heavy-ion reactions at SPS. 
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Fig. 2. Relative differences between the single inclusive pion spectra measured 
in p+p collisions at *fs = 16.9 - 19.4 GeV [11,36,42,43] and the p+p tt +X 
parametrization proposed in ref. [34]. The shaded band represents the 25% overall 
uncertainty that we assign to the parametrization. 



3 Nuclear modification factors at 



20 GeV revisited 



Figure 3 shows the nuclear modification factor, Eq. (1), for 0-7% most cen- 
tral Pb+Pb collisions at ^s NN = 17.3 GeV, obtained using the three p+p 
parametrizations discussed in the previous section. Particle production in 
the low p T region (below p T 1.5 GeV/c) naturally falls below the N coU 
scaling expectation (Raa = 1) since the assumption of independent point- 
like scattering does not hold for soft processes in nucleus-nucleus reactions 
(which instead scale with the number of participant nucleons in the reaction: 
Npart oc [44]). The original WA98 parametrization (open circles) results 
in a steeply rising Pb+Pb over p+p ratio in the whole pt range. The Cronin 
enhancement is apparent above px ~ 2 GeV/c going above Raa ~ 8 at the 
highest Pt's, which is a factor of ~4 larger than the maximum pion enhance- 
ments found in p+A collisions at fixed-target energies (R p a ~ 2) [10,11,12]. 
Wang's p+p reference (crosses) produces an overall less pronounced Cronin 
effect reaching a maximum Raa ~ 3 at high p?. The Raa obtained with the 
Blattnig et al. fit (closed circles) shows no indication of enhancement below 
p T 4 GeV/c and actually, within uncertainties, the data seem to follow the 
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iViV collision scaling (Raa = 1) expected for hard scattering in the absence 
of medium effects. 
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Fig. 3. Comparison of the nuclear modification factors, Eq. (1), for 7r° produc- 
tion in 0-7% most central Pb+Pb reactions at CERN-SPS [7] obtained using the 
three different p+p — ► tt +X parametrizations [7,8,34] discussed in the text. The 
shaded band around Raa = 1 represents the ~10% normalization uncertainty of 
the Glauber calculation of N co u common to all Raa- The bars centered on Raa = 

1 represent the additional fractional uncertainties of the Blattnig (25%) and WA98 
parametrizations (20%) respectively. 

Before addressing the interpretation of the observed high p T pattern in A+A 
collisions at SPS energies, it is legitimate to question the validity of the as- 
sumption of particle production via parton-parton scatterings for such moder- 
ate values of transverse momenta (pt ~ 2.0 - 4.5 GeV/c). Several studies [45] 
have emphasized the relative importance of "soft" effects above p T ~ 2 GeV/c 
in heavy-ion collisions at SPS. Nonetheless, there are at least four indepen- 
dent pieces of experimental evidence which seems to favor an interpretation 
for high p T production in nuclear collisions around y/s NN = 20 GeV based 
on hard scattering processes: (i) All the measured hadron pt spectra above 

2 GeV/c [7,32,33] show a power-law tail characteristic of elementary parton- 
parton interactions, (ii) The shape and width of the near-side azimuthal cor- 
relations of pions with p? > 1.2 GeV/c are jet-like [46], in agreement with 
parton fragmentation expectations, (hi) The measured ratio 7]/tt 0.5 above 
Pt = 1.5 GeV/c [48] is consistent with standard parton fragmentation func- 
tions, and (iv) The observed direct photon yield at px > 1-5 GeV/c [47] is 
consistent with perturbative production cross-sections. 

Using Blattnig's parametrization [34] as proton-proton reference, and the as- 
sociated mean number of collisions (N co u) for each centrality, we present in 
Fig. 4 the nuclear modification factors for high px pion produced near midra- 
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Fig. 4. Nuclear modification factors for pion production at CERN-SPS in cen- 
tral Pb+Pb [7], Pb+Au [32], and S+Au [33] reactions, obtained using the p+p 
parametrization of ref. [34]. The shaded band around Raa = 1 represents the over- 
all fractional uncertainty (including in quadrature the 25% uncertainty of the p+p 
reference and the 10% error of the Glauber calculation of N co u\ (N co u) = 726 ± 
72, 774 ± 77, and 174 ± 20 for Pb+Pb, Pb+Au and S+Au resp.). CERES data 
have an additional overall uncertainty of ±15% not shown in the plot [32]. The 
"curved" band is a theoretical calculation from Vitev and Gyulassy [50] including 
standard nuclear effects (Cronin and shadowing) and final-state parton energy loss 
in a system with initial gluon densities dN 9 /dy = 400 - 600. 



pidity in the three A+A reactions studied at SPS (Table 4) at comparable 
centrality bins. For the three systems, pion production in the range pt ~ 2 - 
4 GeV/c is consistent with Raa — 1- Although with large experimental un- 
certainties, there is an indication of enhancement (Raa > 1) at the highest 
Pt values for the Pb+Pb reaction. The magnitude and pt dependence of the 
nuclear modification factors are compared in the same plot to a theoretical cal- 
culation by Vitev and Gyulassy [50] (yellow band) which includes "standard" 
nuclear effects like Cronin broadening and (anti) shadowing, plus final-state 
partonic energy loss in an expanding system with initial gluon rapidity densi- 
ties dN 9 1 dy = 400 - 600. The influence of nuclear modifications ( "shadowing" ) 
of the parton distributions functions (PDFs) in this kinematical range is small. 
For ^/s NN ~ 20 and pr ~ 2 - 5 GeV/c at mid-rapidity, the colliding partons 
have fractional momenta of the order xsj = 2 Pj ^Js NN = 2 zpx/ y/s NN ~ 0.2 - 
0.4, using (z) = Pt/ P ~ 0.8 for the average fraction of the parton momentum 
P carried by the outgoing (leading) tt°. Such an intermediate x region is mainly 
dominated by valence u, d quarks which are barely modified in the nucleus. 
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The EKS98 parametrization [51] of nuclear PDFs used by Vitev&Gyulassy 
predicts a very modest ~5% antishadowing effect in this p? range. On the 
other hand, the Cronin effect does play a relevant role in 7r° production. The 
Pt broadening is modeled in [50] via multiple scattering in the cold nucleus so 
as to reproduce the magnitude, pt and y/s dependence of the enhancements 
in pion observed in p+A collisions between y/s ^20-40 GeV [10,11,12]. 
For SPS energies the expected effect is R p a ~ 1.4 at pt — 2 GeV/c, steadily 
increasing up to R p a ~ 3 at px = 4 GeV/c (see yellow band in Fig. 5 and 
discussion below). In contrast to this expectation, the central A+A reactions 
at SPS show Raa ~ 1 within p T = 2 - 4 GeV/c (Fig. 4), a factor of ~2 below 
the expected Cronin enhancement. Introducing non-Abelian energy loss [25] 
of the hard scattered partons in a dense expanding system with initial gluon 
densities dN 9 /dy = 400 - 600, provides the suppression needed to reproduce 
the nuclear modification factor observed in the A+A data in this centrality 
bin. Qualitatively similar conclusions have been also reached by closely-related 
pQCD-based calculations of Levai and collaborators [38]. 

Aside from model predictions, the presence of a final-state quenching medium 
in central A+A is confirmed by comparing high-pr tt° production in different 
centrality bins. Fig. 5 shows the nuclear modification factor Raa versus pr 
for peripheral (48%-66% of crpbPb, open triangles), central (0-7% of apbPb, 
circles), and 0-1% most central (closed triangles) Pb+Pb collisions measured 
by WA98. Pions produced in peripheral collisions above px ~ 1.5 GeV/c are 
indeed enhanced compared to "collision scaling", in agreement with the phe- 
nomenological parametrization of the Cronin effect (yellow band) implemented 
in [50]. However, as aforementioned, 0-7% most central reactions are consis- 
tent with Raa being unity up to pr ~ 3.5 GeV/c, and the top 1% most 
central Pb+Pb reactions are actually found to be suppressed, Raa ~ 0.6 ± 
0.15 (stat) ± 0.15 (syst), in this pt range 4 . These results clearly indicate, 
regardless of the p+p — > n°+X reference used in the denominator of Raa, 
that hard hadron production in head-on Pb+Pb collisions at CERN-SPS is 
actually quenched by a factor of ~2 compared to peripheral collisions. It is 
worth to note that this was actually an observation already reported in the 
original WA98 work [7] which, however, remained somehow eclipsed by the 
(conflicting) enhanced values of Raa for all centralities quoted in the same 
papers. 



4 Note that impact parameter fluctuations for this very narrow 0-1% centrality bin 
could result in uncertainties in Raa larger than the quoted ±25% systematic error 
coming from the p+p reference and from the Glauber MC determination of the 
average N coU . 
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19.4 
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Table 1 



Measurements of single inclusive pion production at high pt in heavy-ion reactions 
at CERN-SPS. For each reaction we quote the center-of-mass energy and rapid- 
ity, the experimental rapidity coverage, and the estimated Bjorken energy density 
attained in the most central (0-2%) collisions. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

p T (GeV/c) 

Fig. 5. Nuclear modification factor for ir° production in peripheral (48-66%, (N co u) 
= 78 ± 12, open triangles), central (0-7%, (N co u) = 726 db 72, circles) and most 
central (0-1%, (N co u) = 807 db 81, closed triangles) Pb+Pb reactions at ^/s NN = 
17.3 GeV [7] obtained using the p+p baseline spectrum of ref. [34]. The shaded 
band centered on Raa = 1 represents the overall fractional uncertainty from the 
p+p reference and Glauber calculation of N co u. The curved band is a pQCD-based 
theoretical calculation from Vitev and Gyulassy [50] of pion production in central 
Pb+Pb at 17.3 GeV including standard nuclear effects (Cronin enhancement and 
shadowing) but no final-state parton energy loss. 

4 Discussion 

Figure 6 shows the nuclear modification factors for 7r° production in nucleus- 
nucleus reactions at four different center-of-mass energies and different cen- 
tralizes. Single inclusive pion spectra above pr ~ 2 GeV/c produced at midra- 
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pidity in heavy-ion reactions at SPS and RHIC are suppressed by as much as a 
factor of ~1.6 ± 0.6 (in top 0-1% central Pb+Pb at SPS) and of ~5 ± 1 (in 0- 
10% most central Au+Au at RHIC) respectively, compared to proton-proton 
reactions scaled by the corresponding number of NN collisions. In contrast, 
high pt 7t° production in minimum-bias light-ion (a + a) reactions at ISR 
energies is enhanced (Raa ~ 1-5) with respect to this scaling. Such observa- 
tions are consistent with the expectations of final-state energy loss of the hard 
scattered partons in dense strongly interacting matter produced at midrapid- 
ity in central reactions with heavy nuclei. Determining whether the quenching 
medium is of partonic or hadronic nature (or both) is the ultimate goal behind 
the study of high pt hadroproduction in high-energy heavy-ion collisions. The 
SPS and RHIC measurements of the total transverse energy at central rapidi- 
ties yield values of the Bjorken energy density at r = 1 fm/c, ^ s (t ) ~ 3 
GeV/fm 3 and e™ IC (r ) « 5 GeV/fm 3 , in the top 2% central Pb+Pb,Au [27,28] 
and top 5% central Au+Au [29] reactions resp., which are well above any pos- 
sible scenario involving hadronic degrees of freedom. The equivalent parton 
densities are p SPS (ro) ~ 4.7 fm -3 and p RHIC (ro) ~ 6.9 fm -3 using the thermo- 
dynamical relation p w 2.1 e 3 / 4 (e in GeV/fm 3 ) for a gas of partons as given 
from lattice QCD thermodynamics 5 [1]. The corresponding parton densities 
per unit rapidity, dN/dy — p ■ r • A±, where Aj_ ss 150 fm 2 is the transverse 
area in a head-on A+A collision with A ps 200, are (dN / dy) SPS (r ) ps 700 and 
(dN / dy) nm0 (tq) ps 1100 respectively. Both values are consistent with the ini- 
tial dN 9 jdy gluon densities obtained from the respective "tomographic" parton 
energy loss studies [50] described before 6 . 

Since jet quenching models predict that the induced parton energy loss is 
proportional to the initial parton densities, one could use the measured sup- 
pression at RHIC together with the relative values of esj, or dN/dy(ro), be- 
tween RHIC and SPS to determine the expected value of quenching factor 
Raa at SPS. Indeed, the fact that high-p T hadron spectra in proton-proton 
collisions are well reproduced by a simple power-law l/p? (with n ~ 10 above 
Pt ps 2 GeV/c at SPS, and n ~ 8 above px ~ 4 GeV/c at RHIC), and that 
Raa m central A+A reactions is approximately constant at high-pr both at 
SPS and RHIC, indicates that the basic pr dependence of the yields is not 
changed by the quenching medium (i.e. the power law n exponent remains the 

5 Mind that the use of the ideal gas EOS from lattice calculations with zero bary- 
ochemical potential may be justified at central rapidities for RHIC highest energies 
(where ps << T c ) but is less evident for SPS where ps >T C , and should be taken 
cum grano salis in this latter case. 

6 The values {dN 9 /dy) SPS (T^) w 600 and (dN 9 / dy) Rmc (tv,) 1100 [52] are actually 

/SPS /RHIC 

obtained for r = 0.8 fm/c and r = 0.6 fm/c respectively, for a system with 
transverse area A± ~ 115 fm 2 . However, when scaled to the same area and initial 
proper time (taking into account the decrease in energy density e oc r _4//3 in a 
1+1D Bjorken expansion) one gets coincident results between dN 9 /dyfry) (from 
tomographic studies) and dN/dy(ro = 1 fm/c) (from Bjorken energy densities). 
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8 10 
p T (GeV/c) 

Fig. 6. Nuclear modification factor, Raa(pt), f° r 71-0 production in ion-ion reac- 
tions at CERN-SPS [7] (triangles), CERN-ISR [13] (stars), and BNL-RHIC [15,17] 
(squares and circles). The boxes around the CERN-SPS data points represent the 
normalization uncertainty from the p+p reference and Glauber calculation of N co u. 



same). We can thus relate the suppression in A+A compared to p+p as due 
to a corresponding (energy loss) shift Ap? in the single inclusive spectrum 
dN/dpT- The ratio of A+A over p+p invariant dN/p T dpx spectra, the nu- 
clear modification factor, is then R AA = (1 + Ap T /p T )~( n_1 \ from which the 
corresponding fractional energy loss can be derived: Apx/pr = R AA ~ n ^ ~~ 1- 
Thus, from the measured ~ 0.2 one gets (Apr/p^Wic ~ 0.25, and since 

(dN/dy) SPS « 0.68 (dN/dy) Rmc implies (Ap T ) SPS « 0.68 (Ap T ) RHIC , one would 
expect R S AA 0.4. The minimum value of R A P A measured is, however, R 3 AA xs 
0.6 for the 1% most central Pb+Pb reactions (Fig. 5) which is ~70% larger 
than this simple estimate. The reason for this apparent inconsistency is the 
implicit assumption made above that the counteracting effect of the (Cronin) 
Pt broadening is the same at SPS and RHIC energies. Since the high-p^ spec- 
tra at SPS (with power law exponent n ~ 10) are much steeper than at RHIC 
(n 8) the effect of the initial-state multiple parton scatterings leads to a 
much larger px enhancement at y/s NN ~ 20 GeV than at y/s NN = 200 GeV. 
This fact, as pointed out by Gyulassy and Levai in [30], explains partially why 
the observed suppression at SPS is apparently much lower than at RHIC even 
though the estimated energy densities are only a factor of ~2 larger at RHIC. 

At both center-of-mass energies, however, it is conceivable that not all of 
the high-p^ hadron suppression in central heavy-ion reactions is due to the 
attenuating effects of a partonic medium alone. Indeed, since the produced 
medium undergoes a longitudinal expansion, its initial energy density will 
decrease with time as e oc (r /r) a where a — 1, 4/3 for free streaming and 
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1-D Bjorken expansion 7 respectively [8] in the schematic time-scale evolution 
outlined here. Therefore, even starting at r = 1 fm/c with energy densities 
well above e crit ~ 0.7 GeV/fm 3 , the bulk partonic system will drop below e crit 
and hadronize into a hadron gas phase at r c & r ^f ps 4 fm/c and t^ c ~ 7 fm/c 
respectively. The subsequent hadron system remains strongly self-interacting 
until its density is too low for further rescatterings to take place. At both 
center-of-mass energies, the total lifetime of the strongly interacting system 
(tj ps 15 fm/c [53,54]) is comparable to the time it takes a hard scattered 
parton with typical momentum P = 4 (8) GeV/c to hadronize into a fully 
formed meson [55]: r h ps P ■ R\ ps 12 (25) fm/c (for a pion of radius R h ~ 
0.8 fm [56]). In this context, the produced high-energy parton travels (and 
loses energy) first through a dense partonic system during r < r cr j t and then, 
for a time r crit < r < Tf a , through a hadronic environment. In this second 
hadronic stage, inelastic scattering of the "pre-hadron" object with comover 
soft hadrons of the type described in [5,57] can also partially account for the 
suppression of the final high-pr inclusive spectra. In any case, it is reasonable 
to admit that the energy loss will be larger in the denser (partonic) phase 
than in the more rarefied hadronic one. 



5 Conclusions 

We have reexamined high-pr (pt <^ 2 GeV/c) inclusive pion production in 
nucleus-nucleus reactions from CERN-SPS fixed-target experiments at center 
of mass energies around y/s NN = 20 GeV, and systematically compared them 
to the available proton-proton data in the same range of collision energies per 
nucleon-nucleon pair. In contrast to what has been usually considered so far, 
we conclude that there is no indication of a strong px broadening (Cronin 
enhancement) in the high-j9T yields measured in central Pb+Pb, Pb+Au and 
S+Au reactions. Instead, the data appear to be consistent within errors with 
the perturbative expectations of scaling with the number of nucleon-nucleon 
collisions. The peripheral yields, however, are still found to be Cronin en- 
hanced. These facts, together with the complementary observation that high- 
ly pion yields in head-on (0-1% most central) Pb+Pb reactions are sup- 
pressed by a factor of ~1.6 ± 0.6 compared to p+p collisions, is consistent 
with a moderate amount of final-state quenching of the hard scattered par- 
tons traversing the dense system produced in the course of the most central 
heavy-ion reactions at SPS. Theoretical calculations of parton energy loss in 
an expanding deconfined medium require initial gluon rapidity densities of 
the order dN g /dy ps 600 to reproduce the observed yields, consistent with 

7 The assumption of pure Bjorken (longitudinal boost invariant) expansion is an 
idealistic approximation that can be only applied in heavy-ion collisions at best in 
a narrow rapidity window around y = 0. 



14 



estimations of the Bjorken energy densities attained in the reactions. Addi- 
tionally, we have provided arguments based on hadronization time estimates 
that support the idea that the hard scattered partons must lose their energy 
in the dense partonic and hadronic phases of the reaction. The 2004 Js NN = 




62.4 GeV Au+Au run at RHIC will undoubtedly help to clarify the "excita- 
tion function" evolution of the high-pr hadron suppression observed in high 
energy nucleus-nucleus reactions. The discussion presented here highlights the 
importance of a concurrent and precise measurement of the high-p T produc- 
tion yields in baseline p+p collisions at the same center-of-mass energies as 
the nucleus-nucleus data. 

Note added: Preliminary RHIC results from Au+Au reactions at ^s NN = 62.4 
GeV, obtained after submission of this letter, indicate that high px hadron 
production per NN collision is also significantly reduced (by up to a factor 
of ~3) in central Au+Au compared to p+p collisions at these energies. These 
results are consistent with the indications of moderate high pt suppression in 
heavy-ion reactions at lower SPS energies, discussed here. 
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